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FOREWORD 


This Research Note reports the findings of The llydroloqic Fngineerinn 
Center on appropriate methodoloqies for calculatinq and routinq floods re¬ 
sulting from suddenly-breached dams. 

This study was prepared for the (J.S. Army Engineer Waterways Experiment 
Station, Vickshuro, Miss, with fundino provided bv the Defense Nuclear 
Agency under subtask L10HAXSX337, "Above Ground Structures," work unit 07 , 


"Damage of Dams," and by the Office, Chief of Engineers under DA project 
4A762710AT40, task A1, work unit 006. 

The material contained herein is offered for information purposes onlv 
and should not be construed as Corps of Enoineers policy or as beinn re¬ 
commended guidance for field offices of the Corps of Engineers. 
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RFSFARCH NOTE NO. 5 

GUIDELINES FOR CALCULATING AND ROUTING A DAM-BREAK FLOOn 


1. Introduction . Planning and design requirements for a wide ranqe of pro¬ 
jects, such as emergency preparedness and siting of nuclear power plants, 
have generated widespread interest in dam break floods. Much academic re¬ 
search and some laboratory research have been accomplished on this topic. 
Generalized analytic techniques for calculating and routina such floods, 
particularly in non-prismatic valleys, have not been readily available. 
Furthermore, prototype verification data are almost non-existent. This re¬ 
port describes procedures necessary to calculate and route a dam break flood 
usinn an existing generalized unsteady open channel flow model. The recent 
Teton Dam event was reconstituted to test the model's performance on such a 
highly dynamic wave. The procedures outlined herein relate, primarily, to 
partial breaches. Some deficiencies in the model were identified which will 
require some further research and programming to improve the applicability 
of the program to dam break flood events. 

2 . Summary . The special projects nemo cited as reference (a) established 
four objectives for this study. The first two, a) level of accuracy of 
existing techniaucs and b) sensitivity of calculated results to n-values and 
breach size, are summarized below and presented in detail in Appendix A. The 
third objective, c) description of physical phenomena controllino depth and 
travel time and a discussion of pertinent field data, is presented in the 
body of this report. The fourth objective, d) documentation of the method¬ 
ology, is included in Appendix B. Computer programs utilized in the meth¬ 
odology, references (b) and (c), na.v he obtained from The Hydrologic 
Engineering Center. 

The computer program of reference (c) was applied to the Teton Dam data set 
to demonstrate the level of accuracy one mioht expect in such analyses. 

The results are shown on pages A-?fi through A-?8 of Appendix A and, in 
general, appear reasonable. This test case demonstrates the usefulness of 
a generalized computer program because the methods proposed in references 
(d) and (e) were not applicable to the Teton data set for reasons given in 
paragraph 10. 









Regarding sensitivity to breach size, pages A-22 and A-23 show the two 
breach sizes considered. The breach that developed at Teton was estimated, 
by others, to be 40 percent of the dam embankment. Geometric data were not 
available to verify this, therefore, our best estimate of the final Teton 
breach geometry, page A-22, is based on photographs. The breach shown on 
page A-23 has the same side slope as that on page A-22, 0.6 on 1 but It 
has zero bottom width. This seemed a likely intermediateconditfon, but 
no field data were available at the time of this study to establish an 
observed intermediate condition. 

The calculated outflow is shown on page A-24. The hydrograph labeled 
"trapezoidal breach" assumed the 40 percent breach size, page A-22, 
developed instantaneously. The hydrograph labeled "triangular breach" 
was determined in a similar manner for the TO*? breach size. The third 
hydrograph on page A-?4 was calculated for the trapezoidal breach (labeled 
40?, breach size on page A-' 7 2), but an observed reservoir drawdown curve 
at the dam, page A-20, was used which implies a gradual development of 
the breach rather than instantaneous failure. The last approach was con¬ 
sidered best in estimating the discharge hydrographfrom Teton reservoir 
given the data set and analytical techniques available to us. 

The sensitivity of calculated outflows to breach size and rate of development 
is illustrated on page A-24. It is summarized in the following table together 
with pertinent elevation data for an n value of 0.04. 


Table 1: Sensitivity to Breach Size and Rate of Development 


Final Breach Size 
% of Total Dam 

Rate of 
Development 

[calculated Peak 
Hater Discharge 
at Dam Axis 

106 CFS (1) 

Calculated Peak Elevations 

MSU 

At Dam Axis 



& 1 

171 


T’ r ) 

ITH (3) 

ittuyH 



TO* 

instanta- 

2.4 

5151 

5015 

4933 


neous 





40? 

instanta- 

3.4 

5175 

5020 

4935 


neous 






(1) Multiply by 0.0?032 to get Cubic deters Per Second 

(2) Actual rate of development was unknown so the observed reservoir 
drawdown curve, page A-20, was used to approximate outflow 
conditions. 


















(3) The actual peak discharge, as estimated bv personnel of the Walla 
Jalla District, Armv Corps of Engineers from observed data in the 
Teton Canyon three miles downstream from the dam, was 2,300,000 rfs. 


From these results it is apparent that neither the size of breachs tested 
nor the rates of failure assumed were very significant in predicting peak 
elevations five miles downstream from the dam. 

The calculated peak flood elevations, near the dam, were very sensitive 
to n-values. Increasing n from .03 to .06 raised the peak flood elevation 
25 feet at the dam, as Illustrated on page A-8. T abV At 5 miles 
downstream the calculated difference was only 8 feet. Differences continued 
to diminish with distance. 

Calculated Travel Times are shown on page A-29. They correspond to the 
discharge hydroqraph labeled "simulated from observed data" on paqe A-24 
and n-values of 0.04. 

Searching for a simolified approach in place of references (d) and (e) led to 
a trial application of the Modified Puls routing technique. The hydroqraph 
labeled "simulated from observed data" on page A-24 was routed and a water 
surface profile calculated for the resulting peak discharges. A comparison 
of the results with the observed elevations and the peak elevations computed 
with the full equations is shown on pages A-30 through A-32. Additional 
investigation is neede d to establish the range of applicability of tins * 
method. 


3. Physical Phenomena and Field Data . Analysis of the dam-break flood 
involves understanding the physical processes before aoplyinq analytical 
techniques which approximate those physical processes. Three distinctly 
different processes are involved: the process of structural failure 
causing the breach to develop; the process of setting water into motion 
in a reservoir; and the process of flood wave attenuation. 

The size , shape and rate of breach development are primarily responsible 
for the peak rate of outflow from the reservoir. Yet, of the three physical 
processes, this one is the most difficult to quantify. With the exception 
of man-made breaches, it is difficult to visualize the instantaneous develop¬ 
ment of a breach. Some have occurred, however. The St. Frances Dam, a 
high head concrete gravity structure, apparently suffered an abutment 
failure which resulted in virtually the instantaneous failure of the entire 
structure. The Johnstown flood of 1889 was caused by the complete failure 
of an earth fill dam. Reports Indicate that less than half an hour was 
required for overtopping flow to breach the structure. The recent Teton 
failure, a full depth-partial width breach of an earth fill dam, is estimated 
to have developed in less than two hours. Since natural failure of a major 
structure Is so Improbable, establishing a mode of failure requires a policy 
decision rather than an analytical technique. In qeneral, instantaneous 
failure of the entire structure produces the largest flood wave. 
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The second physical process results from the depth of water above the breach 
invert. That Is, a reservoir has a total enerqy head equal to the elevation 
of the water surface. If the dam Is breached, the force of qravity will set 
water into motion. The effect will propogate, as a negative wave, to the 
upstream end of the reservoir at a velocity equal to yfqy where g Is 
acceleration of gravity and y vs water depth. Because of the great depth 
in a reservoir, very little frictional resistance is mobilized during the 
passage of this negative wave. As a result, water gains specific energy 
rapidly as It moves toward the breach. In instantaneous breach development, 
the peak outflow will occur within a minute or two after breaching . 

Whereas the total energy head setting the water into motion is the specific 
energy (i.e., the initial water depth) above the breach invert, the energy 
which must be dissipated in the downstream channel is equal to the specific 
energy from the downstream channel invert to the initial pool elevation. 

The fact that the water surface elevation drops down rapidly at the dam 
axis does not reflect a corresponding loss in energy head. When flow 
begins, that specific energy above the breach invert is transformed into 
three components: a pressure head, a kinetic energy head and an inertia 
head. (The relative size of each of these energy head components is 
discussed more fully in sections 4 and 5.) Friction loss is relatively 
small and may be neglected unless the reservoir bottom is extremely rough 
(more than 5" or 10" of the water depth). 

The third physical process, flood wave attenuation, involves energy 
dissipation and valley storage. As the flood wave moves downstream, the 
peak discharge tends to decrease, the base of the flood wave will become 
longer and the wave velocity will decrease. Hear the dam, enerqy dissipation 
is primarily responsible for behavior of the flood wave. However, valley 
storage soon becomes the primary factor in flood wave attenuation. The 
key to the transition from energy dissipation to valley storage control is 
the rate at which the slope of the total energy gradient, a line which must 
intersect the initial pool elevation at the dam, is reduced to that of a 
major rainfall flood in the downstream valley. It seems obvious that the 
total energy at any cross section in the valley should not exceed the initial 
reservoir elevation, and yet some analytical techniques occasionally violate 
that constraint. It is good policy to always check the total energy, as 
well as the water volume, in a calculated flood wave. 

The rate of energy dissipation is qoverned primarily by friction loss. 

M inor losses from bends and contractions-expansions are often included in 
the n-values. 

The volume of water in the reservoir is the final piece of field data required. 
This volume strongly Influences the peak elevations at downstream points. 








4. Energy Components and Peak Outflow from Complete, Instantaneous Breaches . 
It is useful to develop the relative size for each energy component in the 
flow at the dam axis and to compare all of them to the more conmon case of 
steady state critical flow at a contraction. 

By assuming a rectangular cross section, zero bottom slope and instantaneous 
removal of the entire dam, Saint-Venant developed an analytical solution for 
the elevation of the free surface, reference (f). page 755. Utilizing that 
equation, the depth of flow at the dam axis was determined, by Saint-Venant 

and others, to be ^ Y 0 where Y 0 is the original water depth at the dam. _ 

Also, the velocity corresponding to the peak outflow was shown to be \ VqY 0 . 
Combining these relationships leads to the equation for peak discharge 


q 


8 

max ’ Tf 



(1) 


Y q is the initial water depth at the dam 
q is acceleration of gravity 
g is peak water discharge in cfs/ft 

m3 X 


Since this equation was developed for a rectanqular section, the total 

discharge may be calculated by multiplying q by the width. 

m3 x 

Using the relationships referenced above, the velocity head (i.e., the 
kinetic energy head component of the specific energy head) was calculated 
to be y- Y„. Since, in the absence of friction and other losses, inertia 
is the only remaining term in the basic, unsteady flow equations of Saint- 
Venant, it may be calculated as follows. 


V ah + — Y + —Y 
o i Q T o g o 


h i 


3 y 
T Y o 


( 2 ) 


These components are shown in Figure 1 along with the energy components 
for critical, steady state flow. 

This figure shows that in the dam break flood analysis, as well as steady 
state critical flow at a contraction, the velocity head is half the pressure 
head. However, the inertia head comnonent is zero in Figure la because flow 
is steady state. 


r. 









(a)Critical flow at a contraction (h)Critical flow from a breached dam 

(steady flow conditions) Unsteady *lo , -< c'w<it1nnsl 

Figure 1. Components of Specific Energy Head. 


The drawdown in water surface elevation to 4" Y 0 at the dam axis, Fiqure lb, 
does not reflect a corresponding energy loss. Experimental results obtained 
by Schoklitsch, reproduced on page 755 of reference If, show relatively 
little friction loss in flow approachinq the dam axis. As might be expected, 
the model results showed friction to be very significant downstream. Tests 
reported hy WES in reference (g)showed no Impact from friction loss at the 
dam axis. However, the WES flume sloped at 0.005 ft/ft, whereas the flume 
in Schoklitsch 1 s experiment had zero bottom slope. 

The significance of this point is that all three enerqy components, pressure 
head, kinetic enerqy head and inertia head, are significant in complete, 
instantaneous breachinqs. Consequently, investiqators encouraqe the use of 
the complete routing equations, often referred to as the Saint-Venant 
equations. Simplifications of the complete equations, such as Musklngham, 
Tatum, Straddle-Staqqer and 'V>dified Puls, are not recommended because the 
empirical coefficients would invariably be developed from rainfall floods 
and would reflect different values of energy components relative to Y 0 . 

5. Instantaneous. Partial Breaches . Partial breaches are classified, 
according to hydraulic performance, as full depth-partial width, partial 
depth-full width or partial depth-partial width. A separate equation has 
been developed for calculating the peak water discharge for each class, 
page 25 of reference (g). 









(3) 


Full Depth-Partial Width 1_ 

0 » JL • h • Y • (£) -v/gY 

max 27 o fr vy o 

B is width of channel, feet 
b is width of breach, feet 


Partial Depth-Full Width j 

'W = 27 • 6 • ^r 1 Yw < a > 

y is depth of water above bottom of breach 


Partial Depth-Partial Width 


'max 


Tf - 


<!> 


1 

Y 1 

f-°) Vi7* 

v v ' 


(S) 


An empirical equation for partial depth-partial width breaches was reported 
in references (g) and (h). 


0 

max 


= 0.29 • h 


y 




• Viy 


( 6 ) 


For breach sizes in the following range. 


1 


<F • r>: ?n 


(7) 


Since the discharge equations for partial breaches are similar, in form, 
to that for a full breach (1), the total specific energy has the same three 
basic components. However, their size, relative to initial water depth, 
is considerably different from that shown in Figure 1. There is no analytical 
solution for partial breaches, therefore, experimental results, presented in 
reference(g), were used to calculate the individual energy head components. 

The following table presents experimental results for full depth breaches 
ranqing in width from 10 % to 100% of the flume width in columns 1, 2 and 3. 
Fractions of initial water depth, calculated with equation 3, are shown in 
columns 4 and 5. A sample of the calculations is presented in the paragraph 
following the table. This sample calculation utilizes equation 3 and a 100% 
breach size (i.e., full breach) to demonstrate that the relative value of each 
energy component is the same as the respective value produced by equation 1, 
the analytical, full breach equation, when equation 3 Is carried to Its 
upper limit. 







Table ?: Relative Size of Enerqy Components in Partial Width Breaches^ 


Test 

No 

Breach 

Size 

(V 

Pressure 

Head 
< Of Y q 

Velocity 

Head (2) 

S of Y 0 

Inertia 
Head (3) 
% of Y 

0 

(1) 

(2) 

(3) 

(4) 

(5) 

1.1 

r ull 

44 

?2 

34 

2.1 

60 

70 

12 

18 

3.1 

30 

8 ? 

12 

6 

4.1 

15 

89 

(4), 

- - 

5.1 

10 

94 

(4) 

- - 


Notes: 1. Values in columns 1 and 2 are from Table A, paqe 8, reference 
(a)and values in column 3 are from experimental results from 
Tables 1 through 5, Station 200, reference (p). 

2. Velocity head is calculated with equations 8 and 9, following. 

3. Inertia head is Y 0 - (pressure head + velocity head). 

4. Calculated values exceeded 100 percent of Y 0 , which probably 
reflects scatter in experimental results. 


wh^re: 

v 

b 


0 

y - max 
max h • y 


= depth of water at dam axis 


= breach width 


0 from equation (3) 

illaX 



For the full breach, b = 1.0B and y = ^ y 

o 


( 8 ) 


(8) 





















1 


„ . ?7 • (of) 1 V?*; 

max 4 

' Y o 

(in) 

* I Viv 

(ID 

v 2 

max = 4 qYo 

2g OTif 

(12) 

" 2 Y 

0 T o 

(13) 


This agrees with section 4 and shows the procedure followed In completing 
Table 2. The inertia head, column 5 in Table 2, was calculated assuming 
zero energy loss upstream from the dam. 


h.- 


i + ? 


+ ? \ 


(14) 



(15) 


Recause of the decrease in relative significance of inertia head and even 
velocity head, it is satisfactory to apply simplifications of the full 
Saint-Venant equations to oartially breached dams. 

6 . Attenuation of the Flood Have . As a flood wave moves downstream, 
friction and other losses chanqe the relative size of the three energy 
components. Even floods from fully breached dams eventually take on the 
characteristics of i rainfall flood and may be routed with a simplified 
routing method such as Modified Puls. Major areas of uncertainty are 1) 
how much distance is required for this transition, 2) how does this distance 
vary when considering partial breaches and 3) what is the maximum breach 
size to consider as a partial breach. 

7. Proposed Analytical Technique . The guidelines presented in Appendix 
R of this report are developed for the computer program "Gradually Varied 
Unsteady Flow Profiles". It is a solution of the basic Saint-Venant 
equations for unsteady flow and may be used to calculate the outflow 
hydrograph through any size or shape of breach, as well as to route that 
hydroqraph downstream and provide water discharge and water surface elevation 
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hydroqraphs at any number of computation points up to 45. The maximum 
discharge, maximum elevation and maximum flow velocity are summarized for 
each computation point. 

Sufficient information is printed out so the time of arrival, time of peak 
and duration of the flood may be plotted. 

This computer program accounts for the movement of the negative wave through 
the reservoir, for the tailwater submergence at the dam, for the three 
components of energy presented earlier, for friction loss and for storaqe 
in the reservoir and the downstream valley. 

Cross sections need not be rectangular or Drismatic. A companion program, 
"Geometric Elements from Cross Section Coordinates", is available to 
transjVm complex cross sections into the required qeometric data set for 
the routing proqram. 

These computer programs are generalized. That is, they are sufficiently 
flexible and adaptable to be used without code changes. They are portable 
from one computer to another and documentation is available, from The 
Hvdroloaic Engineering Center. 

8 . Program Limitations . 

a. Routing with the Gradually Varied Unsteady Flow Profiles computer 
program requires a large high speed computer (50,000, fiO-blt words) and 
personnel who are experienced in applying mathematical models. 

b. Any breach size nay be modeled, but the proqram assumes 
instantaneous development. 

c. All channels must be wet initially. That is, computations cannot 
be made if any portion of the model is dry. This is overcome by prescribing 
a base flow; however, the computer program has difficulty in establishing 
this profile. 

d. ’’ovement of the negative wave through the reservoir causes no 
computational problem until it reaches the upstream end of the reservoir. 
Computation nodes tend to go dry and abort the computer run. 

e. The analysis of multiple failures would require manual intervention 
to stop and restart the calculation process as each new structure is brought 
into the system. 

f. The program assumes a horizontal water surface transverse to the 
flow, whereas a great deal of transverse slope can exist in the actual 
prototype situation. 






y. Pr opose u Areas o f Resear ch. All of the program limitations were 
circumvented in analyzing the Teton Data Set. The trade-off, however, was 
analysis time. Seven weeks were required to set up the data, debug it and 
perform the analysis. The two tasks requiring the most time, probably 75?', 
were establishing initial base flow conditions for the model (8c) and 
stabilizing the computations when the negative surge reached the upstream 
boundary (Bd). Both of these problem areas can be overcome by additional 
programming. The improvements would reduce analysis time to four or 
five weeks. 


Instantaneous breach development, 8b, could be replaced by equations which 
let progressive development take place. In the absence of a theory, the rate 
of development would have to be prescribed with input data. 


Developing the capability to handle multiple dam failures (3e), especially 
in tandem, will he a major modification. 


This analytical technique is a one-dimensional model and will always have a 
Horizontal water surface transverse to the flow. At present, two-dimensional 
modeling is not feasible. 

ID. A1ternate Analytical Proce dures. Alternate analytical procedures were 
proposed in references (d) and Je)■ Neither were applicable to the Teton 
Data Set. 

The dimensionless curves were developed from numerical solution of the St. 
Venant equations and include special treatment of the wave front as it moves 
along a dry channel. By knowing reservoir volume, valley cross section at 
the dam, initial reservoir elevation, stream slope and stream roughness, 
the curves will pro vice three properties of the flood wave: 

1. Time of arrival at downstream points 

2. Maximum depth profile in the downstream channel 

3. Time of maximum depth at downstream points. 

The curves extend for distances ranging up to fifteen times the reservoir 
length. The outflow hydrograph at the dam is not needed to use these 
curves. It was assumed, in developing the curves, that the entire dam 

is breached instantaneously and that the valley is prismatic. Neither 
conoition was satisfied by the Teton case. 

Tne procedure in reference (e)was developed for smaller structures and the 
Teton Data Set was completely beyond the range of nomographs and curves 
presented there. In any case, the procedure does not route the flood wave 
downstream. Only the outflow discharge hydrograph is calculated at the 
uan axis. The procedure can handle a wide range of breach sizes, but it 
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is designed with partial breaches in mind. It has the advantage of tail 
water correction, which is essential when breaching of a low dam coincides 
with a high flow condition In the stream. The procedure is well documented 
and Is easily applied. 

A possible alternative approach for partial breaches is the Modified Puls 
routing technique. Preliminary work with this technique produced the re¬ 
sults shown on pages A-30 through A-32 for the Teton Data Set. A Manning 
n value of 0.04 was used; further details are given In Paragraph 5, Appendix 
A. The advantage of this technique is that readily available and easily 
applied computer programs (e.g., HEC-1 and HEC-2) can be utilized; total 
analysis time would probably be reduced to two to three weeks. 


The disadvantage is that the range of application is limited whereas the 
technique presented in Paragraph 7 is generally applicable. 

Additional research is needed to define the range of applicability of the 
Modified Puls technique. The present hypothesis is that the size of the 
inertia component. Table 2, would provide a suitable parameter for defining 
that range. 

This research would not require additional physical modeling. Studies re¬ 
ported in references (g) and (h) offer test data for numerical studies. 

Other numerical experiments could be performed by using results from anal¬ 
yzing variations of the Teton Data Set with the complete equations. These 
results could be obtained while pursuing any of the areas of research pro¬ 
posed in Section 9. 

Computer programs which utilize the Modified Puls routing technique are 
available and are presently developed to a higher degree of serviceability 
than programs solving the full equations. Water surface profile computations 
will be required In conjunction with the Modified Puls routing to produce a 
water surface profile. These computations are computerized also. No major 
computer program development would be required. The appropriate existing 
computer programs, HEC-1 and HEC-2, are well documented. 
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appendix a 


UNSTEADY FLOW ANALYSIS 
TETON DAM FAILURE 


1. INTRODUCTION 

This phase of the study which calculates and routes the flood resulting 
from a suddenly breached dam, consists of a computer solution in conjunction 
with the Teton Dam failure as defined in objectives a and b of Special 
Projects Memo No. 473. The analysis utilizes the unsteady flow computer 
program to determine water surface elevations resulting from various breach 
sizes and n values. The level of accuracy was determined by comparing avail¬ 
able flood data (particularly high water marks) from the June 5, 1976 dam 
failure with calculated results. 

The analysis used data generally available to field personnel 
engaged in the study of the impact of a dam break flood such as topo¬ 
graphic maps, aerial photography, dam description, gaged stream flows 
and reservoir elevation-capacity curves. The primary area of study 
included about 30 miles of flood plain downstream of the dam. 

2. GENERAL DESCRIPTION 

The Teton Dam is located on the Teton River in southeastern Idaho 
approximately 13 miles north and east of the city of Rexburg (see Plate 1). 
The dam was designed as a zoned earthfill embankment with a crest elevation 
of 5,332 feet (mean sea level datum) and a maximum height of 305 feet 
(above riverbed). It would create a reservoir of 288,250 acre-feet when 
filled to an elevation of 5,320 feet. (Plan and sections of the dam are 
shown in Plates 2 and 3.) The dam is located in a narrow steepwalled 
canyon. The channel geometry is essentially the same for some 4 miles 
downstream of the dam whereupon the Teton River enters upon a wide 
relatively flat flood plain. 

The total reservoir storage just prior to failure on June 5 was about 
251,300 acre-feet at an elevation of 5,301.5 feet. Measured inflow was 
3,580 cubic feet per second and measured outflow before any leaks developed 
was 940 cubic feet per second. 

Although water stored upstream of the dam caused the flooding, other 
sources contributed significantly to the floodflow in the downstream reaches 
However, the complexity of major irrigation diversions and numerous return 
flows precludes an accurate accounting of all sources. According to the 
record obtained from the damaged recorder at the gaging station on Henrys 
Fork near Rexburg ( 






were flowing in Henrys Fork at the time the leading edge of the flood 
wave arrived about 4:00 p.m. on June 5. The contribution of the Snake 
River, just upstream of its confluence with Henrys Fork was estimated 
at 5,600 cubic feet per second around 1:00 p.m. of the same day. 

Table 1 presents a preliminary tabulation of available data pertaining 
to the leading edge of the flood wave generated by the dam failure. 

Distances shown are measured from the dam breach along the estimated 
path of the leading edge of the wave. 

Instantaneous peak discharges were determined at two sites along the 
study reach (Table 2). Indirect measurements based on field surveys and 
empirical formulas were used to compute the peak discharges at the sites. 

The indirect methods for computing peak discharge were based on hydraulic 
equations relating discharge to the water-surface profile and the geometry 
of the channel. 

3. GEOMETRIC MODELS 

The channel configuration upstream and downstream of the dam was 
modeled independently according to procedures contained in the users 
manual, generalized computer program (723-G2-L745B), Geometric Elements 
from Cross Section Coordinates (GEDA), Hydrologic Engineering tenter, (J.S. 

Army Corps'of Engineers, June 1976. 

The upstream geometric model (reservoir) was developed essentially from 
two cross sections located at the damsite and at the estimated upstream 
limit of the reservoir (when filled to capacity). The cross section data at 
the damsite were obtained for preproject conditions from topographic informa¬ 
tion shown on Plates 2 and 3. Information for the other section was obtained 
from U.S. Geological Survey topographic maps (scale 1:24,000). Data were 
coded in the standard HEC-2 1 format and input into the GEDA program. Computed 
results were compared with respective values taken from the reservoir capacity 
curve of Plate 2 and, where necessary, adjustments made to the upstream cross 
section until computed volumes (at specified elevations) plotted reasonably 
close to the published curve (see Plate 4). 

The downstream model (channel) was developed in a similar manner although 
a more detailed definition of the channel configuration was necessary. Since 
the actual flooded area resultinq from the dam failure included several 
major tributaries (the Snake River, Henrys Fork and the Teton River), there 
was no readily available information on storage (volume) within the study 
reach with which to calibrate the geometric model. Cross sections were 
taken on the average every mile, and conveyance limits were designated 
so as to effectively model conveyance for the large flows expected in the 
dam break analysis. 


* HEC-2 Water Surface Profiles , Generalized Computer Program (723-02A), 
Hydrologic Engineering Center, U.S. Army Corps of Engineers, October 1973. 
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TABLE 1 


ESTIMATED TIMES 

FOR LEADING EDGE OF FLOOD WAVE 
(BASED ON FIELD OBSERVATIONS) 


Location 

Miles 3 
downstream 
from dam 

Date 

(month/day) 

Arrival 

time 

(hours) 

Elapsed time 
between sites 
(minutes) 

Approximate 
mean velocity 
between sites 
(ft/min)(mi/hr) 

At damsite 


6-5 

1157 

— 

— 

— 

In Teton Canyon 

3.0 

6-5 


8 

1 ,980 

23 

At Teton 

8.8 

6-5 


25 

1 ,220 

14 

At Sugar City 

12.3 

6-5 


30 

620 

7 

At Rexburg^ 

15.3 

6-5 

1340 

40 


5 

Henrys Fork 
near U.S.G.S. 
gaging station, 
Henrys Fork near 
Rexburg 

22.6 

6-5 

1530 

110 

350 

4 

Menan Bridge on 
the Snake River 
immediately 
upstream of the 
Union Pacific 

30.6 

6-5 

1800 

150 

280 

3 


distances measured along the estimated path of the leading edge of the flood wave. 
b Maximum depths estimated at 6-8 feet within the city of Rexburg. 
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TABLE 2 


ESTIMATED PEAK DISCHARGES 
(BASED ON FIELD OBSERVATIONS) 


Location 

Miles 3 
downstream 
from dam 

Date 

(month/day) 

Time 

(hours) 

Discharge 

(cubic feet per second) 

In Teton Canyon 

3.0 

6-5 

b 


At Teton 

8.8 

6-5 

c 



distances measured along the estimated path of the leading edge of the flood wave, 
kpeak probably occurred between 1230 and 1330 hours. 
c Peak probably occurred between 1300 and 1400 hours. 
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The area, hydraulic radius, top width and average n value were 
calculated at each cross section in the goemetric models for the eleva¬ 
tions specified. By calculating length-weighted values, the preceding 
elements were modified such that they would apply to uniformly spaced 
computation points along the length of each model (reservoir and channel). 

4. UNSTEADY FLOW MODELS 

The unsteady flow data models of the reservoir and the channel were 
developed according to procedures contained in the users manual, generalized 
computer program (723-G2-L2450), Gradually Varied Unsteady Flow Profiles , 
Hydrologic Engineering Center, U.S. Army Corps of Engineers, January 1976. 

Reservoir Model 


In the case of the reservoir, tables of geometric elements were prepared 
at 21 uniformly spaced computation points (approximately 2,800 feet apart) 
and input into the unsteady flow program. Since inflow into the reservoir 
just prior to failure was measured (approximately 3,600 cubic feet per 
second), the upstream boundary condition was based on such. A normal 
depth was determined for the measured discharge rate at the upstream cross 
section. A steady flow condition was assumed, and a constant elevation 
(determined from the normal depth calculation) was specified as the boundary 
condition. Initial values (discharge and elevation) were specified only 
at the most upstream and most downstream computation points (nodes) and 
a linear interpolation scheme within the program utilized to determine 
values at the remaining nodes. The values of discharge and elevation at 
the upstream node are the same as those utilized in the computation of the 
elevation hydrograph at the upstream boundary. The elevation used at the 
downstream node (damsite) was 5,301.5 feet (the observed reservoir water 
surface elevation at time of failure) and the discharge was estimated 
as 3,000 cubic feet per second. (Although the measured discharge was 940 
cubic feet per second before any leaks occurred, a higher discharge rate 
was used because of significant seepage just prior to failure). In the 
reconstruction of the actual discharge hydrograph at the damsite, the 
observed reservoir water surface elevations at particular points in time 
(after failure) were input initially for the downstream boundary condi¬ 
tion (see Plate 5 for a plot of the observed data). Since the observed 
data were not adequate with respect to time to define a reasonable eleva¬ 
tion hydrograph, a greater number of coordinate points were coded in the 
final analysis (refer to Plate 5). The computed discharge hydrograph at 
the damsite, based on the preceding conditions, is shown in Plate 6. (Note 
that the computed maximum discharge is about 1,800,000 cubic feet per 
second at about 2:00 p.m. of June 5. 

When analyzing various breach sizes (partial failures), an elevation 
versus discharge curve was developed for the downstream boundary condition 
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(at the damsite). For the breach sizes shown in Plates 7 and 8 the discharge 
corresponding to various critical depths was calculated according to 1 


where 


Q = 


(b + z D C ) J 
(b + 2z D c ) 



b = bottom width (in feet) 

z = slope of the sides, horizontal divided by vertical 
D c = depth (in feet) 

g = acceleration due to gravity (in feet per second per second) 
Q = discharge rate (in cubic feet per second) 


After converting the various depths to elevations, the discharge rating 
curve corresponding to each breach size was input into the unsteady flow 
program. The computed discharge hydrographs are shown in Plate 9 and 
compared with the simulated hydrograph for observed conditions. 

Preliminary information received indicated that 40 percent of the dam 
embankment was lost such that an initial estimate of the breach configura¬ 
tion (based on available photographs) is shown in Plate 7. The resulting 
discharge hydrograph, Plate 9, represents an instantaneous failure of the 
size depicted. (The same is true in the case of the triangular breach shown 
in Plate 8.) In the case of the discharge hydrograph computed for observed 
conditions, the downstream boundary condition corresponds to a gradual failure 
and results in significant differences in peak and in the time to peak. 

In the unsteady flow analysis of the reservoir system, only the results 
for n values of 0.04 are included. When n values were varied, there was 
no significant difference in the computed discharge hydrographs at the 
damsite. 

Channel Model 


In the case of the channel, tables of geometric elements were prepared 
at 37 uniformly spaced computation points along the study reach (approxi¬ 
mately 4,500 feet apart) and input into the unsteady flow program. The 


handbook of Hydraulics , King and Brater, Fifth Ed., McGraw-Hill 
Book Co., 1963, pp. 8-11. 
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upstream boundary conditions (at the damsite) are the various discharge 
hydrographs computed in the unsteady flow analysis of the reservoir. 

The downstream boundary conditions are discharge rating curves based on 
normal depth calculations. The depth and corresponding discharge was 
determined using Manning's equation 

Q - -±!L_ A R 2 ' 3 S 1 ' 2 
n 


where 

n = coefficient of roughness 
A = area (in square feet) 

R = hydraulic radius (in feet) 

S = slope (in feet per foot) 

2/3 

The values for the geometric elements, A and R , for specified depths 
(elevations), were taken from GEDA output for the most downstream channel 
cross section. The slope was determined from available topographic maps. 
Values of n ranged from 0.03 to 0.07, and for each n value a discharge 
rating curve was developed as a downstream boundary condition. 

In order to simulate the actual flow existing in the probable flooded 
area downstream of the dam, the local inflow option of the unsteady flow 
program was utilized. Major inflow occurred where Henrys Fork and the 
Snake River entered the system. The estimated flows just prior to failure 
were 3,500 and 5,600 cubic feet per second respectively. 

In order to establish stable base flow conditions for each n value 
prescribed, the model was run for a period of 24 hours with the estimated 
inflows at the upstream boundary and local inflow points which existed just 
prior to failure. (The discharge and initial approximations of water 
surface elevation at each node or computation point were input prior to 
running. Since the discharge just prior to failure was estimated from 
observed conditions within the study reach, values were fixed, whereas 
initial water surface elevations were manipulated until stable conditions 
resulted.) Once stable base flow conditions were obtained, the various 
boundary conditions (upstream and downstream) were input into the model. 

The model was then allowed to run until the peak of the flood hydrograph 
passed the downstream limit of study. 

Results of the various runs are tabulated in Tables 3, 4, 5 and 6. 

Plate 10 is a Xerox reduction of the flooded area resulting from the June 5 
dam failure as published by the U.S. Geological Survey. The relative 
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TABLE 4 


COMPUTED AND OBSERVED TIMES 3 
FOR LEADING EDGE OF FLOOD WAVE 


Location 

Miles b 
downstream 
from dam 

Observed 

Arrival 

Time 

(hours) 

0.03 

Computed arrival 
n-values 
0.04 0.05 

(hours) 

■3519 

0.07 

At damsite 

0.0 

— 

— 

— 

— 

— 

— 

In Teton Canyon 

3.0 

1205 

1206 

1207 

1208 

1209 

1210 

At Teton 

8.7 

1230 

1230 

1230 

1230 

1240 

1240 

At Sugar City 

12.9 

1300 

1340 

1400 

1420 

1430 

1440 

At Rexburg 

15.8 

1340 

1440 

1500 

1520 

1530 

1540 

Henrys Fork 
near U.S.G.S. 
gaging station, 
Henrys Fork near 
Rexburg 

21.1 

1530 

1540 

1550 

1650 

1710 

1740 

Menan Bridge on 
the Snake River 
immediately 
upstream of the 
Union Pacific 

27.6 

1800 

1730 

1750 

1830 

1920 

2030 


a Times occurred on June 5, 1976. 

b Di$tances based on the unsteady flow data model. 
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RESULTS OF UNSTEADY FLOW ANALYSIS 
USING SIMULATED DISCHARGE HYDROGRAPH 
FOR 40 PERCENT BREACH SIZE 
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location of the odd numbered computation points and conveyance limits 
are shown (a discussion of the conveyance limits is contained in the 
users manual for the GEDA program). (It is apparent from Plate 10 
that the majority of high water mark observations are along the edges 
of the flooded area.) Maximum water surface elevations are plotted in 
Plates 11, 12 and 13. 

The observed high water marks (taken from Plate 10) were plotted in 
Plates 11, 12 and 13 and designated as left or right bank (looking down¬ 
stream). The ground profile shown reflects the minimum elevation in the 
various cross sections used in the geometric model. Horizontal distances 
are measured along the most probable path of concentrated flow and not 
along any particular tributary as shown in Plate 10. The spacing and 
location of the odd numbered nodes are also included in the profile plots. 
Maximum water surface elevations computed in the simulation of observed 
data were plotted for the odd numbered nodes shown. Water surface profiles 
were then drawn for each n value used in the analysis. 

If a number of cross sections are drawn on Plate 10 normal to estimated 
flow lines, particularly in the area of Rexburg where floodwaters tended 
to spread laterally, a significant variation in high water mark elevations 
between the left and right banks of a given cross section is apparent. In 
some instances, the difference is as much as 30 feet. This situation is 
also apparent in Plates 11, 12 and 13. In general, where floodwaters were 
relatively confined between the left and right banks, elevations of high 
water marks on opposite sides are essentially the same. 

A comparison of the observed and computed maximum water surface eleva¬ 
tions (Plates 11, 12 and 13) would seem to indicate the following: 

(1) Floodwaters were probably concentrated in the center of 
the flooded area between the U.S.G.S. gaging station on 
the Teton River and the downstream limits of Sugar City 
(see Plate 10). 

(2) Floodwaters were probably concentrated in the vicinity 
of Rexburg. 

It is difficult to analyze the accuracy of the computed results based on 
observed high water marks. Any detailed analysis would definitely require 
high water marks within the limits of flooding. If such could be obtained 
a greater effort could be made in stipulating conveyance limits and in the 
actual calibration of an unsteady flow model. 

Tables 1 and 4 compare the estimated (based on field observations) and 
computed arrival times of the leading edge of the flood wave, whereas 
Tables 2 and 5 compare estimated (based on field observations) and computed 
maximum discharges at particular locations within the study area. Plate 14 
gives an indication of the time required for passage of the flood wave 
(based on n=0.04). 
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In the unsteady flow model, determination of the arrival time of the 
leading edge of the flood wave depends on a noticeable difference in the 
computed water surface elevation occurring at a given node. Since minor 
fluctuations of the water surface were present in the base flow and since 
the printout interval was 10 minutes after 12:30 p.m., arrival times were 
estimated (probably ± 30 minutes). 

Computed peak discharges (Table 5) differ significantly from those 
estimated from observed data (Table 2). It is reasonable to assume though 
that the discharges of Table 2 are only rough approximations. 

In the development of the unsteady flow data model, there was no effort 
made to locate computation points at those locations specified in Tables 1 
and 2. (Usually, computation points are located such that results are printed 
at desired locations.) 

5. STEADY FLOW MODEL 

The information originally coded in the HEC-2 format and utilized in 
development of the channel geometric model, was used to determine the storage- 
outflow relationships for various reaches downstream of the dam. Multiple 
profiles were run using HEC-2 and the volumes in a given reach, corresponding 
to a specified discharge, computed. (It should be noted that an n-value of 
0.04 was used in the multiple profile computations.) 

The discharge hydrograph (see Plate 6) based on observed data, was routed 
downstream with the storage-outflow relationships previously determined. 

Channel routing was accomplished with HEC-1 1 using modified Puls. The peak 
discharges computed for each reach were subsequently input into the HEC-2 
data deck. 

Before running, the HEC-2 deck was modified to reflect the conveyance 
limits adopted in the unsteady flow analysis (see Plate 10). Encroachment 
limits were set with the X3 card as close as possible to the conveyance 
limits previously designated (this was necessary since GR stations did not 
always correspond to the conveyance limits specified in the geometric data 
model). The resulting water surface profile is shown in Plates 15, 16 and 17 
(superimposed for comparison). 

The difference in the computed water surface profiles between the two 
methods (the steady flow and the unsteady flow analysis) is due, in part, to 
the difference in the computed discharges. Discharges computed using modified 
Puls were lower in the upstream channel reaches and higher in the lower 
reaches. This is attributed to the significant effect of the inertia component 
(in the unsteady model) inmediately downstream of the dam and in the way 
storage is treated in the two methods. 

1 HEC-1 Flood Hydroqraph Package , Generalized Computer Program (723-010), 
Hydrologic Engineering Center, U.S. Army Corps of Engineers, January 1973. 
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APPENDIX B 


Guidelines for Analyzing a Dam Break Flood with the Computer Program 
"Gradually Varied Unsteady Flow Profiles" 


1. Data Requirements 

a. Topographic maps a re required. The area of coveraqe should extend 
from the upper end of the reservoir in question to the most downstream point 
of interest. The contour interval should be less than the level of accuracy 
expected from the study. 

b. Aerial or ground photographs are desirable. These provide hydraulic 
engineers with a basis for estimating hydraulic roughness coefficients. 
Consequently, they are not needed for the reservoir area but for the channel 
and flood plain. 

c. Reservoir capacity curves or tables are required and may be 
calculated from the topographic maps. These show the volume of water stored 
as reservoir depth increases. 

d. Breach size and shape is required and may be estimated. 

e. Reservoir level at breaching is required and may be estimated from 
aerial photography if reservoir records are not available. 

2. Assumptions . 

a. It is customary to assume that the breach is developed instantly 

and completely to the desired size and shape. Otherwise, several Instantaneous 
breaches of different sizes may be postulated and the resulting floods 
calculated. 

b. Critical depth controls at all partial breaches. In view of the 
values presented in Table 2, main body of this report, inertia is not a 
major consideration over a wide range of partial breach sizes. 

c. The streambed and banks are fixed (i.e., do not erode) during 
the event. 

d. Any bridge across the stream fails instantly upon impact of the 
flood wave. The resulting energy loss is negligible. 






3. General Procedure . The general procedure depends somewhat upon the 
type of breach postulated. That is, if the entire structure is removed, 
the final step of the analysis will calculate the outflow hydroqraph from 
the reservoir plus route that outflow hydroqraph to all downstream points 
of interest in a sinqle operation. 

If only a partial breach is postulated, the outflow hydrooraph is calculated 
first and the routinq to downstream points is accomplished in a separate 
operation. The same analytic technique is used for both. 

The obvious advantaqe of the first approach is a savinqs in the time 
required for performing the analysis. A mor*» siqnificant advantaqe is 
that the method includes any submeroence effects due to the rise in 
tailwater elevation at the dam. This submergence effect can be suffi¬ 
ciently large to change the shape of the outflow hydroqraph from the 
reservoir. 

p artial breaches cannot be analyzed using the first approach because the 
existinq computer program will not accommodate a hydraulic jump. On the 
other hand, the rise in tailwater elevation is usually small enough so 
that no submerqence effect develops. Therefore, the second atmroach is 
adequate. 

4. Calculation of Outflow Hydroqraph Using the Partial Breach Approach. 

The details which follow present the partial breach approach. The first 
step in the calculation will determine the outflow hydroqraph, paragraphs 
a-f. Routing the flood hydrograph to downstream points is presented in 
paragraph 8 . 

a. Calibration of Cross Sectional Data to Reservoir Volume . Locate 
cross sections on the topographic map so reservoir”volume may be calculated 
using the average end area technique. Include all major tributary arms in 
the reservoir. Code this data for the computer program "Geometric Elements 
from Cross Section Co-ordinates" (GEDA) and calculate the reservoir volume 
for the range of elevations up to full pool. This result should agree with 
the reservoir capacity curve. If It differs, ( + 5"), cross sections which 
are more representative should be developed. Usually, this means recoding 
some of the tributary arms. 

b. Geometric Model for Reservoir . Select a distance between computation 
nodes that comnlies with criteria in the routing program, "Gradually Varied 
Unsteady Flow Profiles" (USTrLO). Values typically range between 3 -mile 

(408 m) and 5 miles (15,000 m). Execute GEDA to produce the required tables 
of geometric properties and n-values at each computation node. The most 
downstream section should be at. the dam axis and include the breach that 
was postulated. 
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c. Downstream Boundary . Calculate a critical depth ration curv^ 
for the breach size and shape postulated. Each elevation selected will 
establish an area and top width for use in the followinq equation. 

0 = A -ylq 0 
cr cr ’ cr 


where: 

A „ = cross sectional area, square feet 
cr 


D 


cr 


A /R 
cr 


cr 


R „ = too width, feet 
cr 


o = acceleration of qravitv feet/socond/second 


d. Initial Conditions. Establish the initial water surface elevation 
and discharqe at each computation node in the reservoir. It is exoedient to 
always select a horizontal water surface and zero discharqe. Although, one 
may oostulate whatever conditions he desires bv startinq the comouter model 
at the above conditions and simulatinq inflow/outflow records ud to the time 
when breaching is anticipated. 


e. Upstream Poundary . Usually, the initial reservoir inflow is 
assumed to be zero. It may be otherwise if so desired. In either case, 
code thp inflow hydroqraph as described in the USTFLO users manual. 


f. Calculating the Outflow Hydroqraph . Select an interval between 
printouts that is very short (a minute or less) during the early part of 
the hydroqraph (5 minutes or so). Usually, the calculated peak will occur 
during this interval and the time between printouts may be increased to 1° 
or 15 minutes. In any case, obtain a sufficient amount of Printout to 
define the entire outflow hydroqranh shape. 


5. Routi n') th_e_ bam RreaF Flood 

a. Geometric 'icdel fo_r Houtjni Downstream. Locate cross sections to 
lhp downstream points of interest us inn information oained from the reservoir 
model and followinq s‘eps in the users manual. Select the distance between 
computation nodes. (This distance does not have to be the same as that usM 
in the reservoir when routin') with a two step process. However, use the 
sane qeneral guidelines.) Code tbe data and execute CFHA to produce the re¬ 
quired geometric data tables. 
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b. Downstream Boundary . If possible, set the downstream boundary 
a couDle of nodes beyond the ooint of interest. Calculate a rating curve 
by slope-area or by extrapolating from observed values. 

c. Upstream Boundary . The outflow hydrograph calculated in paragraph 
4 becomes the upstream boundary for this step. 

d. Initial Conditions . Establish initial conditions by calculating 
a base flow profile. (?ero or negative depth is not an acceptable initial 
condition . ) It is expedient to prescribe a bank full water discharge, to 
let the computer model stabilize by simulating the steady flow profile for 
that discharge and finally to decrease the inflow, gradually, until the 
desired base flow discharge is reached. Simulate the base flow conditions 
for a sufficiently lonq period of time to establish a steady flow profile. 
Code this profile (water surface elevation, discharge) to form the initial 
conditions for routing the dam break flood. 

e. n outinfi the Flood . The computer orogram referenced in subparagraph 
lc produces both a discharge hydrograph and a water surface elevation hydro¬ 
graph at every computation nod^. A summary table gives maximum and minimum 
values for elevation, discharge and velocity. 
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